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ABSTRACT 

We investigate in detail some popular cosmological models in light of the latest ob- 
servational data, including the Union2.1 supernovae compilation, the baryon acoustic 
oscillation measurements from the WigglcZ Dark Energy Survey, the cosmic microwave 
background information from the WMAP 7-year observations, along with the obser- 
vational Hubble parameter data. Based on the model selection statistics such as the 
Akaike and the Bayesian information criterias, we compare different models to as- 
sess the worth of them. We do not assume a flat universe in the fitting. Our results 
show that the concordance ACDM model remains the best one to explain the data, 
while the DGP model is clearly disfavored by the data. Among these models, those 
whose parameters can reduce themselves to the ACDM model provide good fits to the 
data. These results indicate that for the current data, there is no obvious evidence 
supporting any more complex models over the simplest ACDM model. 

Key words: cosmological parameters - dark energy - cosmology: observations - 
cosmology: theory 



1 INTRODUCTION 

The present accelerating expansion of the universe is a 
great challenge to our fundamental physics and cosmol- 
ogy. This fact was first discovered by Type la supernov a 
(SNIa) surveys (jRiess et al 1 ll99Sl : IPerlmutter et al.lll999h . 
and later confirmed by precise measurement of the Cosmic 
Micr owave Background (CMB) anisotropics (|Spergel et al.l 
as well as the baryon acoustic oscillations (BAO) in 
the Sloan Digital Sky Surv ey (SDSS) luminous galaxy sam- 
ple l|Eisenstein et al.ll2005t h This cosmic acceleration leads 
us to believe that most energy in the universe exists in the 
form of a new ingredient called "dark energy" which has a 
negative pressure. 

Various theoretical models of dark energy have been 
proposed, the simplest being the cosmological constant A 
with constant dark energy density and equation of state 
wde = p/p = —1- This model, the popular ACDM model, 
provides an excellent fit to a wide range of observational 
data so far. Despite its simplicity and success, the ACDM 
model has two problems. One is the so called "fine tun- 
ing" problem, that is, the observed value of A being ex- 
tremely small comparing with particle physics expectations 
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i|Weinberdli989h . The other is the coincidence problem, i.e. 
the present energy density of dark energy Q,a and the present 
matter density Q m are of the same order of magnitude, 
for no physical obvious reasons. Due to these difficulties 
of the cosmological constant, numerous alternative models 
have been proposed to explain the acceleration instead of 
the ACDM model (see ICopeland. Sami fc Tsuiikawal I2OO6I : 



iFrieman. Turner fc Hutererll2008l for recent reviews) . Gener- 
ally speaking, these models can be divided into two groups: 
one is to modify the matter (i.e. the right-hand side of the 
Einstein equation) and the other is to modify the gravity 
(i.e. the left-hand side of the Einstein equation). 

Although most studies show that the ACDM model is 
in good agreement with observational data, dynamical dark 
energy cannot be excluded yet. In order to distinguish be- 
tween different dark energy models from observations, the 
most commonly used method is to constrain dark energy 
equation of state w. Recent studies have already given tight 
constraints on w, e.g. the Supernova Legacy Survey three 
year sample (SNLS3) combining with othe r probes has given 
W = - 1.061 ± 0.068 (|Sullivan et alj|201lf ). It should be no- 
ticed that although these results are consistent with the 
ACDM model, we cannot yet determine whether the den- 
sity of dark energy is actually constant or whether it varies 
with time as suggested by dynamical models. 



2 K. Shi, Y. F. Huang, and T. Lu 



When one proposes a new cosmological model, it is of 
great importance to place constraints on the model param- 
eters. Usually people use a maximum likelihood estimate 
to set constraints on the parameters of the model. If the 
expected distribution of the data is Gaussian (which is ap- 
plicable for most problems in cosmology), we can use the 
familiar x 2 test for parameter estimation - the smaller the 
X 2 , the better the parameters fits the data. 

On the other hand, in the face of so many different 
dark energy models, a natural question is raised: Which 
model is better, or in other words, which one is most fa- 
vored by the current observational data? This is the problem 
of model selection. One may naively apply the x 2 test here, 
but that does not contain the information of the complexity 
(the number of parameters) of different models. That is, x 2 
statistics are good at finding the best-fit parameters in a 
model but are insufficient for deciding whether this model 
itself is the best one. In order to solve this problem, some 
model selection s tatistics have been proposed in the con- 
text of cosmology ijLiddkfeOoilDavis etall2007l ). The most 
commonly used is the information cri teria (IC) incl uding the 
Akaike information criterion (AIC; lAkaikd 1974) and the 
Bayesian information criterion (BIC; Schwarz 19781 V These 
criterions tend to favor models which give a good fit with 
fewer parameters that embody the spirit of Occam's razor: 
"entities must not be multiplied beyond necessity". 

In this paper, we investigate parameter constraints on 
a number of cosmological models by performing a Markov 
Chain Monte Carlo (MCMC) analysis from latest observa- 
tional data. We then apply the model selection statistics to 
compare different models to assess which is preferred or dis- 
favored by the data. We organized our paper as follows. In 
Section 2 we discuss the model comparison statistics used in 
this paper. In Section 3 we describe the observational data 
used in this paper and the method to use them. Section 4 
gives a detailed description of different cosmological models 
to be tested and the constraining results from observations. 
The comparison of different models by using model selection 
statistics are given in Section 5. The last section presents our 
discussions and conclusions. 



2 MODEL SELECTION STATISTICS 

As mentioned in the introduction, we mainly use the infor- 
mation criteria (IC) including AIC and BIC to test different 
m odels. A detai led description of AIC and BIC can be found 
in lLiddld (|2004l ). The AIC is given by 



AIC= -21n£ max + 2fc, 



(1) 



where C ma x is the maximum likelihood, fc is the number of 
parameters. Note that for Gaussian posterior distribution, 
Xmin = — 2 In Cmax- The AIC was derived from information 
theoretic considerations. 
The BIC is defined as 



BIC = -2\nC max +klnN, 



(2) 



where TV is the number of data points used in the fit. The 
BIC is similar to the AIC, but it includes the number of 
data points in its form while the AIC doesn't. Note that 
for any likely data set In TV > 2, thus the BIC imposes a 
stricter penalty against extra parameters than the AIC. But 



the AIC remains useful as it gives an upper limit to the 
number of parameters which should be included. The BIC 
was derived as an approximation to the Bayesian evidence, 
but this approximation is quite crude. 

The preferred model is the one that minimizes the AIC 
and the BIC. However, the absolute values of them is not 
of interest, only the relative value between different models 
makes sense. 



For the AIC. lBurnham fe Anderson! l|2003h featured the 
following "strength of evidence" in the form of AAIC = 
AIC, - AIC nit n ■ 

AAIC Level of Empirical Support For Model i 
0—2 Substantial 
4—7 Considerably Less 

> 10 Essentially None 



For the BIC, iRobert fe Adrian! (|l995l ) featured the fol- 
lowing "strength of evidence", where ABIC = BICi — 
BIC nil n • 

ABIC Evidence against Model i 

0-2 Not Worth More Than A Bare Mention 
2 — 6 Positive 
6-10 Strong 
> 10 Very Strong 

Thus we can first obtain a model which minimizes the 
ICs, and then we can compare the rest models with it using 
the above judgements as a "strength of evidence". 

It should be noticed that the information criteria alone 
can at most indicate that a more complex model is not nec- 
essary to explain the data, since a poor information criterion 
might rise from the fact that the data are too poor to con- 
strain the extra parameters in the model, and this model 
might be preferred with improved data. 

Furthermore, we must be aware of the limitation of us- 
ing these simplified ICs, since they are based on the best- 
fit x 2 - A more in-depth analysis of model selection should 
consider how much parameter space would give the data 
with high probability, as well as the correlations between 
the parameters. The Bayesian evidence is an approach that 
takes this into account which computes the average likeli- 
hood of a model over i ts pri o r parameter r a nges. S e e, e. 
ISaini. Weller fe Bridle! l|2004h : iLiddlel (|2007f ): iTrottal !i200 
for more discussions. However, the Bayesian evidence needs 
to compute a multidimensional integration over the likeli- 
hood and prior which may be rather complicated. In this 
paper, we'd like to use the ICs instead of the Bayesian evi- 
dence to compare different dark energy models. This simpler 
approach is sufficient for our purpose. 

Besides the ICs, we also apply the reduced chi-square 
and the goodness of fit statistics to see how well the model 
fit the data. The reduced chi-square is xV^i where v is the 
degrees of freedom usually given by TV — fc. It describes how 
well a model fits the observational data sets. The goodness 
of fit (GoF) gives the probability of obtaining a larger dis- 
crepancy between the model and the data than that ob- 
served, assuming that the model is correct. It is defined 
as GoF = r(^/2,x 2 /2)/r(z^/2) where F is the incomplete 
gamma function. 
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3 CURRENT OBSERVATIONAL DATA SETS 

In this section, we describe the latest data sets used in this 
paper and the method to analyze them. 



3.1 Type la Supernova (SNIa) 

Currently, SNIa is the most powerful tool to study dark en- 
ergy because of their role as standardizable candles. 
For the SNI a data, we use the currently largest Union2.1 
compilation (|Suzuki et al.ll201lT ) that contains a total of 
580 SNI a, which is an updated version of the Union2 com- 
pilation (lAmanullah et al1l2010l ). The newly added twenty 
supernovae are all at relative high redshift (0.6 < z < 1.4) 
and thus can help tighten the constraints on the evolution 
behavior of dark energy. 

Cosmological constraints from SNIa data are obtained 
through the distance modulus y,{z). The theoretical distance 
modulus is 



Hth(zi) = 51og 10 D L (zi) + no, 



(3) 



where fj,o = 42.38 — 5 log 10 h with h the Hubble constant Ho 
in units of 100 km/s/Mpc, and the Hubble-free luminosity 
distance Dl is defined as 



Vl(z) — — sinn 



(4) 



where E(z) — H(z)/Ho and Q.k is the present curvature 
density. Here the symbol sinn(x) stands for sinh(x) (if fi*. > 
0), sin(x) (Q k < 0) or just x (fl k = 0). 

To compute the y 2 for t h e SN Ia data, we fol- 
low iNesseris &: Perivolaropoulosl |2005f ) to analytically 
marginalize over the nuisance parameter Ho- 



100 h" 1 Mpc separation measured from the Sloan Digital 
Sky Survey (SDSS) Third Data Release (DR3) Luminous 
Red Galaxy (LRG) s ample with effective redshift z = 0.35. 
IPercival et al.l (|2010h performed a power-spectrum analy- 
sis of the SDSS DR7 dataset, considering both the main 
and LRG samples, and measured the BAO signal at both 
z — 0.2 and z = 0.35. Recently, in the low-redshift uni- 
verse the 6dF Galaxy Surve y (6dFGS) team ha s reported 
a BAO detection a t z = .1 (|Beutler et alj|201lf ). Most re- 
cently, iBlake et al.l |201ll ) presented measurements of the 
BAO peak at redshifts z = 0.44,0.6 and 0.73 in the galaxy 
correlation function of the final dataset of the WiggleZ Dark 
Energy Survey. They combined their WiggleZ BAO mea- 
surements with SDSS DR7 and 6dFGS datasets to give tight 
constraints on dark energy. In this work, we follow them to 
constrain different dark energy models using their combined 
BAO dataset. We highlight our usage of this combined 
BAO dataset, since there are altogether six data points, 
which are more than previous BAO data, and few have 
used this combined BAO datasets to constrain dark energy 
since the publication of the WiggleZ paper. 

The data can be found in the above papers, but for 
completeness here we summarize the BAO measurements 
and the way to use them. 

The y 2 for the WiggleZ BAO data is given by 
IBlake et al.l |201ll ). 



XwigglcZ — (A bs ~ Ath)CwigglcZ (A bs ~ Ath) 



(8) 



where the data vector is A obs = (0.474, 0.442, 0.424) for the 
effective redshift z — 0.44, 0.6 and 0.73. The corresponding 
theoretical val ue Ath denotes the ac oustic parameter A(z) 
introduced by lEisenstein et all (120051 ) : 



Xsn = A - 2/io-B + nlc, 



(•>) 



A{z) 



D v {z)^/a^Hl 



(9) 



where 



.4 



B 



= E 

i=l 
580 

= £ 

580 



fj-obs(zi) — nth(zi-, no = o) 



(6) 



C : 



a is the uncertainty in SNIa data. Eq. (5) has a minimum 
for /xo = B/C at 



-2 

Xsn 



A c ■ 



(7) 



This equation is independent of fj-o, so instead of Xsn we 
will adopt Xsn to compute the likelihood. 



3.2 Baryon Acoustic Oscillations (BAO) 

The competition between gravitational force and primordial 
relativistic plasma gives rise to acoustic oscillations which 
leaves its signature in every epoch of the universe. As stan- 
dard rulers, BAOs provide another independent test for con- 
stra ining the property of d ark energy. 

lEisenstein et al.l (|2005l ) first found a peak of this baryon 
acoustic oscillations in the 2-point correlation function at 



and the distance scale Dv is defined as 



Dv(z) = — 
no 



-,1/3 



(10) 



where Da(z) is the Hubble-free angular diameter distance 
which relates to the Hubble-free luminosity distance through 
Da{z) = Dl(z)/(1 + z) 2 . The inverse covariance IeZ 
is given by 



r -i 

^ Wig 



1040.3 -807.5 336.8 \ 
-807.5 3720.3 -1551.9 . (11) 
336.8 -1551.9 2914.9 / 



Similarly, for the SDSS DR7 BAO distance meas ure- 
ments, the x can De expressed as l|Percival et al.ll2010h 

Xsoss — {dobs — dth)C S o SS {dobs — dth) T , (12) 

where d b s = (0.1905,0.1097) is the datapoints at z = 0.2 
and 0.35. d t h denotes the distance ratio 



d z = 



r a (z d ) 

Dv(z) 



Here, r s (z) is the comoving sound horizon, 

C -^ldz' 



r.iz = c 



H(z') 



(13) 



(14) 



where the sound speed c s (z) = l/\/3(l + Rb/(1 + z), with 
R b = 31500n b h 2 (T C MB /2.7K)" 4 and T C mb = 2.726K. 
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The redshift za at t he baryon drag epoch is fitted with 
the formula proposed by lEisenstein fc Hul |l998), 



Zd = 



1291(fi m /i 2 ) c 



1 + 0.659(fi m /i 2 ) - 828 



where 



61 = 0.313(f7 m /i 2 )-°' 419 [l + 0.607(tt m /i 2 ) ' 674 ] 
b 2 = 0.238(f3 m /i 2 ) 0223 . 



(15) 



(16) 



Csdss m Eq. (12) is the inverse covariance matrix for 
the SDSS data set given by 



CI 



30124 -17227 
-17227 86977 



(17) 



For the 6dFGS BAO data (jBeutler et al.ll201lf ). there is 



only one data point at z = 0.106, the x is easy to compute: 

(18) 



2 i'd z - 0.336 

XsdFGS — 



0.015 



The total \ for all the BAO data sets thus can be 
written as 

XBAO = XWiggleZ + XsDSS + XedFGS- (19) 



3.3 Cosmic Microwave Background (CMB) 

Since the SNIa and BAO data contain information about 
the universe at relatively low redshifts, we will include 
the CMB information by using the WMAP 7-yr data 
l|Komatsu et al.ll201ll ) to probe the entire expansion history 
up to the last scattering surface. The \ 2 f° r t ne CMB data 
is constructed as 



Xcmb — X T C c Ij B X, 



where 



X = 



l A - 302.09 
R - 1.725 
z* - 1091.3 



Here Ia is the "acoustic scale" defined as 

7rdx,(z*) 



Ia 



(l + z)r a {z,)' 



(20) 



(21) 



(22) 



where cIl( z) = Dl(z)/Hq a n d the redshift of decoupling z, 
is given bv lHu fc Surivamal l| 19961 ). 



2* = 1048[1 + 0.00124(f7 6 /i 2 )~°' 738 ][l + gx{Q. m ti z ) 



2\921 



9i = 



0.0783(fi b /i 2 



1 + 39.5{n b h 2 ) c 



,92 = 



0.560 



l + 21.1(fi i ,/i 2 ) 1 - 8 



(23) 



(24) 



The "shift parameter" R in Eq. (21) is defined as 
(|Bond. Efstathiou fc Tegmark|[l997h 



R = 



;Dl{z) 



c(l + 0*)~ 

Cq MB in Eq. (20) is the inverse covariance matrix, 



^CMB — 



2.305 29.698 -1.333 
29.698 6825.270 -113.180 
-1.333 -113.180 3.414 



(25) 



(26) 



Table 1. Summary of cosmological models 



Model 


Abbreviation 


Parameters 


Cosmological constant 


ACDM 




Constant w 


roCDM 


fife, Q m , w 


Varying w (CPL) 


CPL 


U k , tt m , Wo, W a 


Generalized Chaplygin Gas 


GCG 


U k , A s , a 


Dvali-Gabadadze-Porrati 


DGP 




Modified Polytropic Cardassian 


MPC 




Interacting Dark Energy 


IDE 


f2fc, f2 m , w x , <5 


Early Dark Energy 


EDE 


Hfe, flm, fl e , WQ 



Note: The Hubble constant Ho in the fit is not deemed as a 
model parameter, but we include it in the number of degrees of 
freedom and in k when calculating the AIC and BIC. 

3.4 Observational Hubble Data (OHD) 

In addition to the SNIa, BAO and CMB data, we also use the 
observational Hubble parameter as an observational tech- 
nique. These data compose an independent dataset that can 
help break the parameter degeneracies, thus may also shed 
light on the cosmological models we aim to study. 

In this work, we adopt 11 data points from differential 
ages of old passive evolving galaxies l|Stern et al.ll201oT ). the 
\ 2 value for these OHD can be expressed as 



2 

XOH 



[H th { Zi ) - H obs { Zi )f 



(27) 



where Oi is the la error in the OHD data with Zi ranging 
from 0.1 to 1.75. 



4 COSMOLOGICAL MODELS AND 
CONSTRAINING RESULTS 

In the following, we study eight popular cosmological mod- 
els discussed in the literature. The models with their pa- 
rameters and the abbreviations we use are listed in Table 
1. We examine them through the expansion history of the 
universe to see whether they are consistent with current 
data at the background level. The model parameters are 
determined through the minimum \ 2 fitting by using the 
Markov Chain Monte Carlo (MCMC) method. Our MCMC 
code is based on the publicly available CosmoMC package 
l|Lewis fc Bridle! 120021 ). 

It should be stressed here that unlike most other 
work on dark energy model constraints, we do not as- 
sume a spatially flat universe as a prior in this paper, al- 
thou gh recent studies sh owed that the universe is nearly 
flat dKomatsu et al.ll201lf) . When we constrain the proper- 
ties of dark energy, the parameters such as the equation of 
state w, are always degenerate with the curvature density 
Qt- It has already been shown that ignoring Qk will in- 
duce large errors on the reconstructed dark energy param- 
eter (e.g. w). If the true geometry is not flat and with the 
wrong flatness assumption, one will erroneously conclude a 
wrong behavi or of dark energy even if the curvature term 
is very small (|Clarkson. Cortes fc Bassetdl2007l: I Zhao et alj 
120071 : IVirev et al.ll200Sft . So in our work, instead of assum- 
ing a flat universe, we will include i} k as a free parameter in 
different cosmological models. 



A comprehensive comparison of cosmological models from latest observational data 5 



4.1 Cosmological constant model 

The cosmologi cal constant A was originally introduced by 
lEinsteirJ (|l917n to achieve a static universe but later aban- 
doned by Einstein after Hubble's discovery of the expansion 
of the universe. Ironically, after 1998 the cosmological con- 
stant revived again as a form of dark energy responsible for 
the late-time acceleration of the universe. The cosmological 
constant plus cold dark matter (CDM) is usually called the 
ACDM model, and in this model the dark energy equation 
of state w = — 1 at all times. The Friedmann equation in 
this case is 

h 2 {z)/hI = o(i+z) 2 +O(i+.z) 4 +0n(i+z) 3 +(i-O-O-fi 

(28) 

where the radiation density parameter O is given by O- = 
n 7 (l+0.2271AT e// ) with fi 7 = 2.469 x 10~ 5 fo ~ 2 and the effec- 
tive n umber of neutrino species N e j f — 3.04 (Ko matsu et alj 
l201ll ). We caution that in many papers the O term is usu- 
ally neglected. While this is reasonable for SNIa analysis 
where the redshift is very small, for high redshift especially 
at CMB epoch this term is dominated. When calculating 
the sound horizon r s for CMB and BAO analysis, ignoring 
this radiation term will induce large errors on the results, so 
it should better be included. The last term in the equation 
represents the energy density of the cosmological constant. 

This simple model has only two parameters O and Q, m . 
Our global fitting from all the four data sets gives the best-fit 
values with 1 a errors: 

O = -0.0024 ± 0.0056, On = 0.291 ± 0.014. (29) 

Our results are co nsistent with the latest results of the 
WiggleZ BAO paper l|Blake et alJ^OHl ). Fig. 1 shows the 
constraint from each of the SNIa, CMB and BAO data sets 
and the joint constraint from all the four data sets. We do 
not separately give the constraint from OHD data since cur- 
rently it is not as stringent as the first three probes, but 
we include it in the combined results. It can be seen that 
although the contour of each single data set is quite broad, 
their combined constraint is quite stringent and this reminds 
us of the power of joint analysis from different independent 
data sets. A flat universe is quite favored by current data 
within la confidence level. 

4.2 Constant w model 

The simplest extension to the ACDM model is to assume 
that the dark energy equation of state w does not precisely 
equals —1, but a constant to be fitted with data. In this 
model, the Friedmann equation is 

H 2 (z)/H 2 = + zf + 0(1 + zf + 0.(1 + zf 

+ (l-O-O-0n)(l + ^) 3(1+tu) . (30) 

There are three parameters in this model: O, f2m, and 
w. The best-fit values using all the data sets are 

O = -0.0012±0.0064, On = 0.292±0.015, w = -0.990±0.041 

(31) 

also in agreement with lBlake et all l^OHl ). 

We plot the contour of Q. m and w after marginalizing 
over O and Ho in Fig. 2. This model also gives a good 
fit to different data sets. The combined result shows a clear 



preference around the cosmological constant model (w = — 1 
within la confidence level). 

4.3 Chevallier-Polarski-Linder model 

There is no prior reason to expect w to be —1 or a con- 
stant, and if w varies with time, the Friedmann equation is 
modified as 

H 2 (z)/Hl = 0(1 + z) 2 + 0(1 + zf + On(l + zf 

+ (1 - O - O - On) exp ^3 ^ 1 ' dz ^j ■ ( 32 ) 

) ? 

Many function forms of w evolving with reds hift have 
been proposed so far (see, e.g. IJohri fe Rath]|2007l ). Among 
the various parametrizati ons of dark energy equa t ion o f state 
w, the one d eveloped by IChevallier fc Polarskil (|200lh and 
iLinderl ((2003) turns out to be an excellent approximation 
to a wide variety of dark energy models, and this CPL 
(Chevalier-Polarski-Linder) model is the most commonly 
used function form to study the time dependence of w. The 
equation of state in this model is 

w(z) = wo + w a ^ ^ , (33) 

So insert Eq. (33) into Eq. (32), we get the Friedmann 
equation for this CPL model: 

H 2 (z)/Hl = 0(1 + zf + 0(1 + zf + Q m (l + zf 

+ (1 - O - O r - On) (1 + Z f (1+W0+W ^ exp(^^). 

1 + z 

(34) 

There are four parameters in this model: ilk, £^ m , wo 
and w a - Our best-fit values for these parameters are 

O = 0.00027lo;oosi, On = 0.293 ± 0.016, 

Wo = -0.966_o.io5, w a = 0.202I 1 5 3 . 

Fig. 3 shows the contours of wo and w a for the CPL 
model after marginalized over other parameters. Obviously, 
w a is weakly constrained by current data. This is partially 
due to the degeneracy between the curvature and the equa- 
tion of state. If we set O = in the fit as most work did, 
the constraints would be more stringent, especially for a 
single data set. However, as explained earlier, we don't as- 
sume a flat prior in the fitting procedure. We see once again 
that it is consistent with the ACDM model f or wp = — 1 
and w a = 0. Our results are in agreement with Blake et alj 
(|201lT ) although they assumed a flat universe in their fit. 

4.4 Modified Polytropic Cardassian expansion 

The Cardassian ex pansion model was first proposed in 
iFreese fc Lewisl l|2002t ) which modifies the Friedmann equa- 
tion to allow for an acceleration in a matter-dominated uni- 
verse. The motivation for this modification could be the em- 
bedding of our observable universe living as a 3-dimensional 
brane in a higher dimensional universe. The original form of 
the Cardassian model can be written as 

H 2 (z) = ^ Pm + BpZ, (36) 

where B is a constant and n is a dimensionless parameter. 
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This power law form is equivalent to the constant w 
model (Sec. 4.2) for w = n — 1, so there is no need to ad- 
ditionally fit this model. Here we consi der a modified poly - 
tropic Cardassian model proposed by Wa ng et al.l {2003), 
and in addition, we also include the curvature and radiation 
term : 



H 2 {z)/H 2 Q 
+ft m (f+z) 3 



Q k (l + z) 2 +Q r (l + z] 

1 — — Si] 



1+ ( 



s2n 



)«-!(! + *) 



3q(n-l) 



(37) 

The above equation reduces to the ACDM one for q = 1 
and n = 0. Our joint constraints give the best-fit parameters 
as follows: 



Sl k = 0.0022 ± 0.0025, tl m = 0.280 ± 0.006, 

n 00-7+0-152 n c «o+0.S56 

q = 0.897_ . 468 , n = -0.648_i.m6. 



(38) 



The constraints on the parameter q and n is very weak 
from current data. Fig. 4 displays the marginalized contours 
of q and n. It can be seen that it is still consistent with the 
ACDM model in la level. 



4.5 Dvali-Gabadadze-Porrati model 

The Dvali-Gabadadze-Porrati (DGP) model is a popular 
model which modifies the grav ity to allow for cosmic accel- 
eratio n without dark energy (|Dvali. G abadadz e fe Porratil 
2000). This model may arise from the brane world theory 
in which gravity leaks out into the bulk at large scales. The 
Friedmann equation is modified as 



h 2 {z)/hI = a k (i + zf + r« r (i + z) 



+ (yn m (i + Zf + n ra + 



(39) 



where r c is the length scale beyond which gravity leaks out 
into the bulk, and Q. rc = l/(4rlH%). Setting z = in Eq. 
(39), we get the normalization condition 



(1 



4(1 - flfc - Q r ) ' (40) 

The DGP model has the same number of parameters as 
the ACDM model. The marginalized best-fit parameters are 
as follows: 



n k = 0.020 ± 0.006, n m = 0.305 ± 0.015. 



(41) 



We can see that although the matter density is con- 
sistent with that of the ACDM one, the curvature term is 
much larg er than that in other mode ls. This feature i s also 
noticed bv lZhu fc Alcanizl 1120051 ) and lGuo etafl (|2006h who 
got a non-fiat universe for the DGP model at high confidence 
level. In Fig. 5, it can be seen that the three observational 
probes strongly disagree - the areas of intersection of any 
pair are distinct from other pairs. CMB data prefer a posi- 
tive J}fc_whiIe_SN and B AO data are in s uppor t of negative 
f2 fe . lRubin et al.l l|2009f) and lDavis et al.l (|2007r i also noticed 
this signal. It may imply that this DGP model is strongly 
disfavored by current data. This can be further quantified 
by the model selection statistics to be shown in Sec. 5. 

It should be mentioned that the DGP model could per- 
form better when using only SNIa data. For example, us- 
ing MLCS2k2 light curve fitter for SDSS-II supernova data, 



ISollerman et"aH l|2009fl found that the DGP model perform 
better than the ACDM model under the information crite- 
ria. Also, recently it was noticed that the Supernova Legacy 
Survey (SNLS) three years data analyzed wi th SALT2 fit- 
ters a lone prefer the DGP model over others (|Li. Wu fc Yul 
I2012F ). However, when combining with BAO and CMB data, 
things have changed, and the concordance ACDM model 
became favored. This is not surprising, since current SNIa 
data are mainly confined by systematic errors rather than 
statistical errors, so it would be better to combine SNIa data 
set with other probes (BAO, CMB, etc.) to constrain cos- 
mological models to avoid any potential bias that may be 
caused by the systematics of SNIa. 

4.6 Interacting Dark Energy model 

The fact that the energy density of dark energy is the same 
order as that of dark matter in the present universe sug- 
gests that there may be some relations between them. This 
may rise from the an interaction between a scalar field (e.g. 
quintessence field) and dark matter. Such motivation may 
help alleviate the coincidence problem. A popular approach 
to study this interaction is to introduce a coupl ing term on 
the r ig ht hand side o f the c ontinuity equations jDalal et al 
200ll: ICai fc Wand liool; iGuo. Ohta fc Tsuiikawal |200 ' 



Caldera-Cabral. Maartens fc Urena-Lo pcz 2009): 



f)ra + 3Hp m = +Tp m , 

p x + 3H(1 + w x )p x = -Yp n 



(42) 



where F characterizes the strength of the interacting, p m is 
the matter density and p x is the dark energy density with 
w x the equation of state. In order to place observational con- 
straints on the coupling term T, it is convenient to express 
F in term of the Hubble parameter H 



T = 5H, 



(43) 



where 5 is a dimensionless coupling term. Note that a posi- 
tive 5 corresponds to a transfer of energy from dark energy 
to dark matter, whereas for a negative 8 the energy transfer 
is opposite. 

It is obvious the expansion history will depend on the 
parameter S , and thus we are interested in placing observa- 
tional constraints on it. For simplicity, here we assume 8 to 
be a constant. In more general case 8 may be varying, and 
there have already been a lot of work on this varying case. In 
this paper since we mainly focus on the model comparison, 
studying a constant coupling is enough for our purpose. 

For a constant 8, solving Eq. (42) with Eq. (43), the 
Friedmann equation becomes 

H 2 {z)/H% = fi fe (l + z) 2 +fir(l + z) 4 + (l-n m -fi fc -n r )(l + 2 



,3(l+ii; je) 



On 



[6(1 + z) 



3(l+iu x ) 



+ 3w x (l + z) 



3-6 



(44) 



8 + 3w x 

This model has four parameters: Sl k , ^m, w x and 8. The 
concordance ACDM model is recovered for 8 = and w x = 
— I. Our global fitting gives the following best-fit values 



fife = 0.0007 ± 0.0032, Q m 
8 = -0.0043 ± 0.0066, w x 



0.292 ± 0.007, 
-1.001 ± 0.087. 



(45) 



Fig. 6 shows the case for this interacting dark energy 
model. It is noticed that the SNIa and BAO data sets give 
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quite weak constraints on the parameter space, comparing 
to the CMB data set. This is not strange since the SNIa and 
BAO data are located in low redshifts, and we can see from 
Eq. (44) that when z <C 1 the 5 + 3w x term just cancels 
out, so the corresponding information about 8 lost. This 
tells us how important it is to include other high redshift 
data. The CMB and OHD data are appropriate for this 
purpose. Our results show that the ACDM model still re- 
mains a good fit to the data (at least within 2 a level), 
but a negative coupling (5 < 0), i.e., the energy transfers 
from dark matter to dark energy, is slightly favored. Also in 
this case the equation of state of dark energy w x prefers a 
phantom case w x < — 1. Thi s result is consistent with that 
obtained by other authors llGuo. Qhta fc Tsuiikawa 20071; 



Valiyiita. Maartens fc Maierottol l20ld ; ICao. Liang fc Zhul 
201ll ). 



4.7 Generalized Chaplygin Gas model 

Despite of their quite different properties in the equation of 
state and clustering, the temptation to unify dark energy 
and dark matter in a single entity has occured to many re- 
searchers from the beginning. To realize it, a natural and 
simple way is to introduce a perfect background fluid. The 
Chaplygin gas model is just a typical example. 

The original Chaplygin gas model w as proposed by 
iKamenshchik. M oschcll a fc Pasquierl ([2001) . In this model, 
the pressure P of the fluid is related to its energy density 
p through P — ~A/p where A is a positive constant. In a 
more general case, one may consider a generalized chaplygi n 
gas model (GCG) given bv lBento. Bertolami fc Sen! (2002) 



-Ap<*. 



(46) 



Consider the energy conservation in the framework of 
Friedmann-Robertson- Walker (FRW) metric, we obtain the 
following solution 



p{a) = po 



A s + 



A s 



a 



3(l + a) 



(47) 



where A s = A/py +a , and po is the present energy density of 
the GCG. One finds the intriguing feature that the energy 
density of this GCG acts like dust matter in the early time 
and behaves as a cosmological constant at late epoch. So 
the GCG model can account for both dark matter and dark 
energy at the background level. The Friedmann equation for 
this model can be written as 

H 2 (z)/Ht = 0(1 + zf + 0(1 + zf + 0(1 + zf 

+ (1 - O- - Or - n b )[A 3 + (1 - As)(l + z) 3(1+a) ]T^, 

(48) 

where O is the present density parameter of the baryonic 
matter. We adopt Q b = 0.0451 ac cording to the WMAP 7- 
year results ()Komatsu et al.ll201ll ). The effective total mat- 
ter density can be expressed as On = O + (1 — O — O — 
0-)(l-A) 1/(1+Q) . Note that the concordance ACDM model 
is recovered by a = 0, thus On = 1 — O — O — A s (l — O — 
Ofc — O). 

There are three parameters in this model. The fitting 
results are: 



fife = 0.0004±0.0032, As = 0.733±0.025, a 



-0.011±0.140. 

(49) 



The GCG model provides a good fit to the data. Fig. 7 
shows the contours for the two parameters ^4,, and a in this 
GCG model. It can be seen that the ACDM model (a = 0) 
falls well within la level, and the original Chaplygin gas 
model (a = 1) is ruled out at more than 2 a confidence level. 
This i s in agreement with the results of iLiang. Xu fc Zhul 
Hon]) and lWu fc Yul l|2007f) . 



4.8 Early Dark Energy scenario 

One of the differences between dynamical dark energy and 
the cosmological constant is that the energy density of the 
former may be non- negligible even at very high redshift (e.g. 
around recombination, or earlier). The existence of the so- 
called "tracker" field jSteinhardt. Wang fc Zlatevlll999h is 
important to alleviate the cosmological constant problem. 
The tracker fields correspond to attractor-like solutions in 
which the field energy density tracks the background fluid 
density for a wide range of initial conditions. These models 
can be motivated by dilatation sym metry in particle physics 
and string theory (Wetterichlll988h . 

As a specific model of such early dark energy scenario, 
here we consider a co mmonly used form with t he dark energy 
density expressed as |Doran fc Robbers! [iool l 



0.de{z) 



O 



0[l-(l + z) 



(1 + z)- 



+ O e [l - (1 + zf 



(50) 

where Q.%, E is the present dark energy density, O is the 
asymptotic early dark energy density and rao is the present 
dark energy equation of state. Thi s equation is based o n sim - 
ple considerations as depicted in Doran fc Robbers! (|2006l ) 
and Doran. Schwi ndt fc Wetterichl l|200ll ). The early dark 
energy behavior is included in the O term. The —Swo term, 
motivated by the relation Ubb(z) /Q m (z) oc (l + z) iw , allows 
the deviation from the ACDM model. 

Eq. (50) assumes a spacially flat universe, and since in 
this work we do not assume flatness from the beginning, 
we would like to slightly modify that equation to include a 
contribution from curvature, 

n DE (z) = (q° de - o[i - (i + ^) 3 ™ ]) / 

[^DE +On(l + z)-' iWa + O (1 + Z)- 3W0 + 1 + 0(1+ Z)-^ - 1 ] 

+ Oe[l- (l + z) aW % (51) 

where fi'p-E = 1 — — O — O- 

In this case, the Friedmann equation can be expressed 

as 



r2/ . wo .2 _ On(l + Z) 3 + O(l + -Z) 4 + O(l+Z) 2 



1 — Q,de{z) 



(52) 



This model also has four parameters. The best-fit values 
are as follows 



Qk = 0.0042 ± 0.0069, On = 0.291 ± 0.007, 
O = 0.026±o.o26> ™o = -1.039 ± 0.097. 



(53) 



As we can see from Fig. 8, the ACDM model (O = 
0, wo = —1) is still favored within 2a level. However, the 
early dark energy component is not totally excluded from 
current data. Since the SNIa and BAO data sets are low- 
redshift ones, they cannot give effective constraints on the 
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early dark energy density fi e , so the most stringent con- 
straint comes from th e CMB data s et. Th e results are con- 
sistent with those of iRubin et al.l l|2009h . ICalabrese et al,l 
l|201lf) and lReichardt et alj l|2011 



5 MODEL COMPARISON 

In this section, we compare the above different models by us- 
ing the model selection statistics. Table 2 gives a summary of 
the IC results. It is easy to see that the concordance ACDM 
model has the lowest ICs, so the AAIC and ABIC are all 
calculated with respect to the ACDM model. 

Given current data sets, the ACDM model is clearly 
preferred by these model selection tests. Following it are 
a series of models that give comparably good fits but 
have more parameters. According to their ICs, we can 
roughly rank these models to four groups: 1. positive against 
(GCG,wCDM), 2. strong against (IDE,EDE), 3. very strong 
against (MPC,CPL), 4. essentially no support (DGP) from 
current data. The GCG and uiCDM model fit the data well, 
maybe due to their fewer parameters and that they can eas- 
ily reduce to the ACDM model. The IDE and EDE models 
are punished by the ICs mainly because they have more 
parameters. The constraints on the MPC and CPL mod- 
els are very weak, and they are also penalized by their large 
number of parameters. We see that the DGP model is so 
strongly disfavored by the data as its AICs have much larger 
values than others. Its goodness of fit is also much smaller 
than others. So we can say, at least at the background level, 
the DGP model can be excluded by current joint data sets 
at high significance from a model selection point of view. 

To see more clearly how to realize cosmic acceleration 
from these models, we plot the deceleration parameter q in 
Fig. 9. The deceleration parameter q, defined as q — —da/a 2 , 
can be calculated by 

1 + 2 dH(z) 



q = -l + 



H{z) dz 



(54) 



As expected, these models all give negative q at late 
times, and positive q at earlier epoch, meaning that the ex- 
pansion of the universe slowed down in the past and speeded 
up recently. Phenomenologically, there is a transition red- 
shift z t between the two epochs, and we also give it in the 
figure. We can see from Fig. 9 that much due to their com- 
plexity, the constraints on the CPL and MPC models are 
very weak, as the contours of their parameters. Although 
the constraint on the DGP model is quite tight, it gives 
the transition redshift z t = 0.45, much smaller than the 
other models, suggesting a strong distinction of the DGP 
model from other models. Given the bad behavior of the 
DGP model from the model selection techniques discussed 
earlier, this smaller transition redshift, too, may suggest that 
the DGP model is disfavored by current data. The concor- 
dance ACDM model remains the best fit in the figure. 



6 DISCUSSIONS AND CONCLUSION 

In this work, we have studied a number of different cos- 
mological models in light of the latest observational data. 
The data we used include the newly published Union2.1 su- 
pernovae compilation and the WiggleZ BAO measurements, 



together with the WMAP 7-year distance priors and the 
observational Hubble data. By using these data sets, we ob- 
tained the best-fit parameters for different models. We use 
the information criteria including the AIC and the BIC, to 
compare different models and to see which is the most fa- 
vored one by current data. These ICs tend to favor models 
that give a good fit with fewer parameters. Unlike many pre- 
vious work did, we do not assume a spatially flat universe 
in our work, instead, we treat the spatial curvature fife as a 
free parameter in the fitting procedure. 

Using the AIC and BIC for model comparison, it is 
found that the concordance ACDM model remains the best 
one to explain current data. The generalized chaplygin gas 
model and the constant w model also give good fits to the 
data. The interacting dark energy model, the early dark en- 
ergy scenario, the Chevallier-Polarski-Linder model and the 
modified polytropic cardassian model are all punished by 
their large number of parameters, thus are not favored by 
the ICs. The DGP model gives the worst fit, although it 
has the same number of parameters as the ACDM model. 
Its AIC and BIC are much larger than other models, with 
a bad goodness-of-fit. Meanwhile, the curvature den- 
sity parameter fife is quite near zero for all models 
except the DGP model. In the DGP model the dif- 
ferent contours from the different observational data 
sets strongly disagree— SNIa and BAO prefer nega- 
tive value of fife whereas the CMB prefers positive 
fife, and so the joint constraint on the value of fife is 
much larger than in other models. 

We also showed the deceleration parameter q for dif- 
ferent models, finding that all models indicate a late time 
cosmic acceleration consistent with observations. However, 
the transition redshift zt for the DGP model is much smaller 
than that in other models. This may reflect the fact that 
the DGP model can not reduce to the concordance ACDM 
model for any value of its parameters. 

In brief, given current data sets, the ACDM model re- 
mains the best one from a model-comparison point of view, 
followed by those that can reduce to it. Those who can not 
reduce to the concordance model fit the data quite badly, 
especially for the DGP model. In spite of its observational 
success, due to the theoretical considerations, we can not 
yet say that this ACDM model truly describes our universe. 
For the time being, we can at most conclude that this model 
fit the current data best among various models. With more 
and more precise data available in the future, it is excepted 
that we will finally be able to identify the nature of cosmic 
acceleration. 
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Figure 1. The marginalized la and 2a contours of the ACDM model parameters and £7 m , derived from different data sets. "ALL" 
denotes the joint constraint including all the four data sets. 
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denotes the joint constraint including all the four data sets. 
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11.83 


EDE 


555.06/595 


87.81 


3.08 


11.87 
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DGP 
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13.01 



Note: The ACDM model is preferred by both the AIC and the 
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Figure 3. The marginalized la and 2a contours of the CPL model parameters wq and w a , derived from different data sets. "ALL" 
denotes the joint constraint including all the four data sets. 



Liddle A. R., 2004, ??jnlMNRAS, 351, L49 
— , 2007, ??jnlMNRAS, 377, L74 
Linder E. V., 2003, ??jnlPbys. Rev. Lett., 90, 091301 
Nesseris S., Perivolaropoulos L., 2005, ??jnlPhys. Rev. D, 
72, 123519 

Percival W. J. et al., 2010, ??jnlMNRAS, 401, 2148 

Perlmutter S. et al, 1999, ??jnlApJ, 517, 565 

Reichardt C. L., de Putter R., Zahn O., Hou Z., 2011, 

ArXiv e-prints 
Riess A. G. et al., 1998, ??jnlAJ, 116, 1009 
Robert E. K., Adrian E. R., 1995, Journal of the American 

Statistical Association, 90, 773 
Rubin D. et al., 2009, ??jnlApJ, 695, 391 
Saini T. D., Weller J., Bridle S. L., 2004, ??jnlMNRAS, 

348, 603 

Schwarz G., 1978, The Annals of Statistics, 6, 461 
Sollerman J. et al., 2009, ??jnlApJ, 703, 1374 
Spergel D. N. et al., 2003, ??jnlApJS, 148, 175 



Steinhardt P. J., Wang L., Zlatev I., 1999, 
??jnlPhys. Rev. D, 59, 123504 

Stern D., Jimenez R., Verde L., Kamionkowski M., Stan- 
ford S. A., 2010, ??jnlJ. Cosmology Astropart. Phys., 2, 

8 

Sullivan M. et al., 2011, ArXiv: 1104. 1444 

Suzuki N. et al, 2011, ArXiv e-prints 

Trotta R., 2007, ??jnlMNRAS, 378, 72 

Valiviita J., Maartens R., Majerotto E., 2010, 

??jnlMNRAS, 402, 2355 
Virey J.-M., Talon-Esmieu D., Ealet A., Taxil P., Tilquin 

A., 2008, ??jnlJ. Cosmology Astropart. Phys., 12, 8 
Wang Y., Freese K., Gondolo P., Lewis M., 2003, ??jnlApJ, 

594, 25 

Weinberg S., 1989, ??jnlRev. Mod. Phys., 61, 1 

Wetterich C, 1988, Nuclear Physics B, 302, 668 

Wu P., Yu H., 2007, ??jnlApJ, 658, 663 

Zhao G.-B., Xia J.-Q., Li H., Tao C, Virey J.-M., Zhu 



12 K. Shi, Y. F. Huang, and T. Lu 




Z.-H., Zhang X., 2007, Physics Letters B, 648, 8 
Zhu Z.-H., Alcaniz J. S., 2005, ??jnlApJ, 620, 7 



A comprehensive comparison of cosmological models from latest observational data 13 




Figure 5. The marginalized la and 2a contours of the DGP model parameters Q m and f2fe, derived from different data sets. "ALL" 
denotes the joint constraint including all the four data sets. 
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Figure 6. The marginalized lc and 2<r contours of the IDE model parameters <5 and w x , derived from different data sets. "ALL" denotes 
the joint constraint including all the four data sets. 
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Figure 7. The marginalized la and 2a contours of the GCG model parameters A s and a, derived from different data sets. "ALL" 
denotes the joint constraint including all the four data sets. 
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Figure 8. The marginalized la and 2a contours of the EDE model parameters wq and Q e , derived from different data sets. "ALL" 
denotes the joint constraint including all the four data sets. 
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